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Abstract 

We describe a new method based on algebraic techniques, which leads to a model 
of atom-diatom collisions. 


1 Introduction 

One-dimensional atom-molecule collisions have been studied using a combination of differential 
and algebraic techniques in conjunction with time evolution operator methods for a variety 
of molecular potentials [1]. For three-dimensional systems, however, one needs to resort to a 
number of computational integro-differential techniques to describe such interactions [2,3]. In 
recent papers we have proposed an algebraic framework based on the vibron model [4], that 
leads to a three-dimensional scheme for such interactions [5], which starts from a U( 4) x U ( 4 ) 
description of three-atomic molecules. We then apply a coherent state method to one of the 
U{ 4)’s, thus extracting a coordinate dependence for the interaction between one of the atoms and 
the remaining diatom, while the latter is still described algebraically. The resulting Hamiltonian 
is formally analogous to the ones used in one dimension [1] and can be solved in principle in the 
interaction picture [5]. This is a difficult task, however, particularly with respect to assessing the 
different degrees of approximation involved. For this reason we consider here a one- dimensional 
version of our model based on a 17(2) x 17(2) dynamical algebra, which is simpler to analyze and 
already incorporates much of the complexity of the full three-dimensional model. 


2 One-Dimensional Model 

In our one dimensional model, the stretching vibrations of triatomic molecules are described in 
terms of the dynamical algebra 


U x (2)xU 2 (2) 


( 1 ) 



by means of the Hamiltonian 


H = E 0 {N\ ,N 2 ) + AJ 2 Zl + BJ] 2 + CJ Z1 J t2 + D J x ■ J 2 + . . . , (2) 

where J Xi , J yi and J Zi are the SU t (2) generators and N, the (fixed) total boson number associ- 
ated to Ui(2). The dots at the end of (2) indicate that we may also need other combinations of 
generators, such as J yi • J y3 . Such Hamiltonians are well suited to describe stretching vibrations, 
due to the connection that can be established between the U (2) algebra and the one-dimensional 
Morse oscillator Schrodinger equation [6]. We now introduce the coherent state basis associated 
to bond number 2 [7], 


| [N]r >= 


1 

v/An(l + r 2 )" 


(s' + rt') N |0> , 


( 3 ) 


where s are two scalar boson creation operators through which the U (2) algebras are realized 
[6]. Computing the expectation value of (2) in the basis (3) and carrying out the coordinate 
transformation 


g — 6r/ao 


r = 


-bx/ao 


( 4 ) 


we arrive at a Hamiltonian of the form 


where 


and 


H = £ 0 + t>, + V 2 + V 3 , (5) 

H 0 = AJ 2 ti , (6a) 

V, = -de _6z/2ao ( 2 - e -6l/ao ) 1/2 , (66) 

V 2 = -P(2e~ bx/ao - e - 2hxo/ao ) , (6c) 

t> 3 = 7e _6l/ao , (6 d) 

a = -±(C + D)NJ Zl , (7a) 

p = -±BN(N- 1) , (76) 

7 = -\DNJ yi . (7c) 


The potentials (6b) and (6c) generalize the typical interactions of exponential form (6d) used 
in previous works [1]. For collision energies high compared with the vibrational excitations, we 
may “freeze” the molecular coordinates and substitute J 2 , and J yi in (7a) and (7c) by their 


expectation values in the | jm > basis, leading to a well defined potential for which we evaluate 
classical trajectories, associated to a classical Lagrangian of the form 


L= fi i 2 +V 1 (x) + V 2 (x) + V 3 (x) . 
£ 


( 8 ) 


Substituting 
the form 


these trajectories x(£) back into (5) then leads to a time-dependent potential of 


V(t) = v, + V 2 (t) + V 3 (t)J yi + V 4 (t)J Xl , 


(9) 


where Vi(t) are obtained from (6), except for V 4 (t), which is identically zero. 

We now use the interaction scheme, which is appropriate to determine time-dependent 
solutions for Hamiltonians of the form H = H 0 + V(t), as in (5). The evolution operator in this 


scheme satisfies the equation 




( 10 ) 


where 

V ini (t) = e ihHo{i - tl) V(t)e- ihHo{t - tl) . (11) 

A great simplification 1 arises when Vi n t is a linear combination of a closed algebra. For the 
particular case of SU( 2), can be expressed in the form 


= C~ e ” l /^53(0'lz g — 


( 12 ) 


If this is the case, substitution of (12) into (10) gives rise to differential equations for gi(t) in 
terms of the (known) potentials Vi(t). Once is known, the scattering matrix is defined 

in terms of it as 

S = U(- 00 , 00 ) , (13) 

i.e., the transition probabilities can be obtained through 

P(\jm >— > | jm >) = | < > | 2 

r 771 2 


where — <jfj( — 00 , 00 ^/% and the d-functions are the usual Jacobi functions appearing in the 
definition of Wigner’s D-functions [8]. 


3 Calculations 

The main stumbling [block for the evaluation of the evolution opperator is the calculation of 
Vi nt (t) through (11). We proceed to discuss briefly some approximations tb carry out this task. 
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a) Sudden Approximation 

In the sudden approximation, the interaction time is considered to be very short. A Taylor 
expansion of (11) along t\ leads to 

VS»*(i) = V(t) + l -[H 0 ,V(t)]V(t -U) + ... (15) 

which at zeroth order gives Vi nt (t ) = V(t). This represents an approximation where all dy- 
namical information of the molecular potential is lost in equation (10). We may then solve for 
the classical trajectories using (12), in some cases analytically, which generalize results known 
previously for the exponential interaction (6d) [9]. In this limit we are able to exactly reproduce 
the analysis of Levine and Wulfman [10,11] and extend them to the more general interactions 
in (6). We omit these results here for lack of space and refer the reader to reference [9]. 

t 

b) Averaging Techniques for Vi n t(t) 


Since we are interested in comparing our algebraic methods with quantum mechanical ab 
initio calculations (when these are available), we should improve the evaluation of V, nt so as to 
introduce the dynamical effects of the molecular potential. The sudden approximation of the 
previous subsection fails to account for the fine properties required. Fortunately, a quantum 
mechanical calculation for the transition probabilities in the scattering of an atom from a one- 
dimensional Morse oscillator has been presented by Clark and Dickinson [12], so we may gauge 
our approximations by comparing with their results. Returning to eq. (11), we substitute (6a) 
and make use of the commutation relations 

[J], Jy] = ~J y ~ 2iJ z J x 

[J z , J z\ — —J X +2 iJ z Jy 


(16) 


After some algebra we find the expression 


V, nt(t) = /" — -r ( — cos(2x 0 b} 0 tj z )jy + J y cos(2 x 0 u> 0 tj z j) 

2 sinfxoc^ot) ' ' 

+ - - — T ~ r (- cos(2x 0 uj 0 tJ z )J y + J y cos(2x 0 u> 0 t J*)) \V 3 (t) 

2sin(xou;of) ' / J 

+ V l (t)J Zt +V 2 (t) , 


(17) 


where xo^o is the standard anharmonicity parameter, which can be expressed in terms of the 
Morse potential parameter through xqojo = The problem with eq. (17) is that it no longer 
satisfies being linear in the SU( 2) generators, so (12) is not valid. To de^l with this problem, 
we proceed as follows. We apply V,„t(f) to a Morse eigenfunction |jm > and insert a complete 
set of states 1 3 1 1 >< jt 1 1 t° the resulting expression. Taking into account the selection rules 


for Jy, which imply, /x = m ± 1, we find a state-dependent potential: 


v m (t) 


-( 


cos(2ar 0 ^o^ m ) 

cos(x 0 ct;o<) 


V 3 (t) + 


cos(2x ou> 0 t(m + 1)) 

cos(xQ<ji>ot) 


+ 1 >< 3 m + 1| 


+ !)> v 3 (t)l jm - 1 X jm - 1|)) j, + v,(t)l + V 2 (t) 

cos(xoo;o<) / 


(18) 
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The m-dependent cosine functions, however, are practically the same for small values of the 
phonon number v = j — m (large m values), so we can simplify (18) by averaging the excitation 
to the two possible states /i = m ± 1, so we arrive at the simple m-dependent potential 

V m (t) = (cos(2xou>otm) cos(x 0 wo*)^3(0^y 

+ Vt(t)J z + V 2 (t) , (19) 

which is linear in the )>£/( 2) generators, so (12) is valid. Since we are interested in comparing 
our approach with the calculations of ref. [12] where only Vj(t) ^ 0, the potential (19) simplifies 
further and only J y remains. The evolution operator is then given by 


U m {t) = 


( 20 ) 


where 


g 2 (m,t) = Vicos(xouot)cos(2xouomt) . (21) 

The differential equation for the evolution operator then lead to exactly solvable forms for the 
classical trajectories and to the asymptotic value of g 2 {m,t) (see below (14)) 


E f°° 2 ( [~E t \ 

/3 2 (m) = -jT J cos(x 0 u> 0 t)cos(2mx Q u>o) sech 


(It 


( 22 ) 


where E is the collision energy. From (14) this leads to the simple expression for the transition 
probabilities 

P(\jm >— > | jm' >) = d } m ,(/3 2 (m)) , (23) 


where 


n ( t\ _ I ~ | v ’ " I 

2 m / ^ 2 ^ [sinh((2m' — l)xoko) sinh((2m' + l)zoko)J 


(2 m' — l).r 0 


( 2 / 77 ^ 1)^0 


(24) 


with k Q = In Table I we compare the results of our calculation, using (23), (24), 

with those of Clark ahd 

TABLE I. Comparison between the algebraic model and the Clark-Dickinson Model 12 * 

Energy ( hu>o / 2 units): 16 


Initial State = 0 

Final state Clark-Dickinson 
Probability 

1 0.245000 

2 3.38000E(-02) 

3 2.72000E(-03) 

4 1.61000E(-04) 


Initial State = 1 


Model 

Probability 

0.157837 

1.45652E(-02) 

8.65153E(-04) 

3.71652E(-05) 


Clark-Dickinson 

Probability 

0.318000 

5.72000E(-02) 

5.16000E(-03) 


Model 

Probability 

0.253196 

3.62444E(-02) 

2.86219E(-03) 
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5 


7.44000E(-06) 1.22993E(-06) 


3.24000E(-04) 1.50963E(-04) 


Energy ( huio/2 units): 12 

Initial State = 0 Initial State = 1 


Final state 

Clark-Dickinson 

Model 

Clark-Dickinson 

Model 


Probability 

Probability 

Probability 

Probability 

1 

0.129000 

0.124019 



2 

6.65000E(-03) 

8.57432E(-03) 

0.167000 

0.208828 

3 

1.49000E(-04) 

3.81572E(-04) 

8.41000E(-03) 

2.19917E(-02) 

4 

1.48000E(-06) 

1.22807E(-05) 

1.38000E(-04) 

1.28983E(-03) 

5 

4.95000E(-09) 

3.04486E(-07) 

7.01000E(-07) 

5.07090E(-05) 


Energy ( hu>o/2 units): 10 


Initial State = 0 Initial State = 1 

1 7.39000E(-02) 0.105782 

2 1.66000E(-03) 6.09133E(-03) 8.15000E(-02) 0.182403 

3 1.17000E(-05) 2.25776E(-04) 1.22000E(-03) 1.58582E(-02) 

4 1.95000E(-08) 6.05219E(-06) 5.12000E(-06) 7.71373E(-04) 

5 0. 1.24982E(-07) 2.95000E(-10) 2.51949E(-05) 

i 

Energy ( huo/2 units): 8 

Initial State = 0 Initial State = 1 

1 3.15000E(-02) 8.66178E(-02) 

2 1.85000E(-04) 3.98810E(-03) 2.28000E(-02) 0.152903 

3 1.18000E(-07) 1.18193E(-04) 3.07000E(-05) 1.05375E(-02) 

4 0. 2.53330E(-06) 1.37000E(-11) 4.08112E(-04) 

5 0. 4.18294E(-08) 0. 1.06318E(-05) 

Energy ( hu>o/ 2 units): 6 

Initial State = 0 Initial State = 1 

1 7.31000E(-03) 6.64920E(-02) 

2 2.30000E(-06) 2.29488E(-03) 1.25000E(-03) 0.120126 

3 0. 5.09821E(-05) 0. 6.15327E(-03) 

4 0. 8.19113E(-07) 0. 1.77900E(-04) 
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5 


0 . 


1.01384E(-08) 


0 . 


3.46550E(-06) 


Dickinson [12] for different collision energies and for two different initial states. Although our 
results differ from the exact ones, particularly for higher final states where the probabilities are 
very small, the general trend is reproduced remarkably well, taking into account our semiclassical 

method. 


4 Conclusions 

The algebraic model seems to provide an attractive alternative to integro-differential techniques 
for the description of atom-molecule collisions. The approximation methods developed for the 
one-dimensional test model can be readily extended to the three dimensional case and applied to 
real systems [2,3]. We are currently exploring the generalization of these techniques to include 
reaction channels, which would represent an important development because of the relevance of 
these processed in atmospheric interactions [13]. 
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